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o : 

^ Dimensional arguments lead to say that excitons in a sample of volume V behave as bosons 

for \Na^j./V <^ 1, Ox being the exciton radius and A a dimensionless factor. The Mott criterion, 
' which is based on the disappearance of all exciton bound states because of screening, gives A ~ 1. 

Since excitons feel each other through both Coulomb interaction and Pauli exclusion between their 
electrons and holes, criteria based on the underlying fermionic character of the exciton should be 
QQ ' even more relevant. Two significant quantities are the boson departure l — [B, B^] and boson number 

B^ B, where B^ is the exciton creation operator. Their expectation values in the A''-exciton state 
are close to their boson values for A ~ 100 and A ~ 50, respectively. By contrast, the expectation 
CJ I value of B^ (B^)^ in the vacuum state is close to iV! for 25 iV^ al/V <.l only. This surprising A^^ 

■ dependence comes from the intrinsic A'^-body character of Pauli exclusion. Consequences of these 

new criteria on the Bose condensation of excitons are discussed. 



When the density of electron-hole (e-h) pairs in semiconductors is very large, the carriers form an e-h plasma. 
Although the plasma's Coulomb energy is large, it remains small compared to its kinetic energy so that, in this 
C regime, the Coulomb interaction can be treated as a perturbation (with possible high-order processes, such as RPA 
I or "bubble" diagrams, in order to avoid spurious divergences The situation is quite different in the low-density 
, regime, where one electron basically interacts with only one hole, giving birth to a hydrogenoid "atom" called exciton. 
O ' Since bound states only appear if the Coulomb interaction between electron and hole is treated exactly, i.e., to all 
^ O^ , orders in perturbation, an exciton may be viewed as the result of a repeated interaction between electron and hole 
through "ladder" diagrams 

^-H ' From dimensional arguments, we expect the cross-over between the two regimes to appear when the distance 
I> between two excitons is of the order of their size, i.e., when XNa^/V ~ 1, Ux being the exciton Bohr radius, TV the 
^^O number of excitons in the sample volume V , and A a dimensionless factor. 

, The high-density regime in which the Coulomb energy is dominated by the kinetic energy is valid when the parameter 
' Ts, defined by 1 7r(rs GxY/V — 1, is small compared to 1. This leads to A = 47r/3 ~ 4. The prefactor A can also be 
obtained by evaluating the carrier density for which, due to Coulomb screening, no excitonic bound state survives. 
' The accepted result |Qj for this Mott density leads to A ~ 1. 

[ According to the spin-statistics theorem, excitons, which are composite particles made of two fermions, should 
"Y^ ' behave as bosons. While the above A's have been deduced from considerations in which the Coulomb interaction 
^ plays a crucial role, criteria based on the underlying fermionic character of the excitons should be even more relevant. 
^ ■ Indeed, excitons are not perfect bosons because they "feel" each other not only through Coulomb interaction but also 
' through Pauli exclusion between their electrons and holes . 

In this Letter, we address the question: Up to what density can excitons be treated as bosons? We show that the 
criteria based on the fermionic character of excitons generate prefactors A larger than 1 by two orders of magnitude. 
Q i This casts new light upon the possible observation of the challenging Bose condensation of excitons. 

' For the sake of simplicity, we will restrict here to excitons with zero total momentum. In terms of the free electron 
(aj.) and hole (&Lj,) creation operators, the exact creation operator for the ground-state exciton reads 



X3 
O 



Bt=^0,45L„ (1) 

k 

where 0^ is the ground-state exciton wave function in momentum space. From eq. (|l]), we get 

D = 1-[B,B^=Y, \ A\Halak + bl,^b.k)- (2) 

k 

Obviously, D would vanish for perfect bosons. It is usually said that this operator, which measures the departure 
from bosonic behavior, "is of the order of N a'^/V" . Such a claim pertains to an operator expectation value only, so 
that one must specify the state on which the expectation value is taken. For the A^-free-pair state 
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\i^k,...k.) = alb\---al^bl,jO), (3) 

we find 

N 

ii^k, ...k.mk, = 2 J2 I'Pkf -2N\q^k=o\' (4) 

if all ki^s are close to zero. For V sufficiently larger than — so as to avoid boundary effects — the bulk exciton 
ground state gives |(/)fc|^ — ^/(l + a^)'*, where z — GAira^/V. The expectation value of D on is thus less 

than 1 for A = 1287r « 400. 

Actually, the relevant state in this problem is not the A^-free-pair state but the iV-pair state constructed from the 
ground-state exciton 

|^,»,, . _ (g')"l»> ,5, 

^{0\B^ (St)^|0} 

This state is the valid picture for the ground state of N pairs in the low-density limit, as shown by Keldysh and 
Kozlov 1^. As the quantity 

{0\B^{B^)^\0)^NIFn (6) 

appears in any expectation value taken on j'l'*^^^), we will start by its calculation. 

If B were a perfect boson operator, we would get = 1. From eq. (p|) we find instead 



F2 = al - (T2 

-F3 = (Ti - 3 CTi 0-2 + 2 0-3 

F4 — af ~ 6 (jf (T2 + 8(Ti (T3 -I- 3 cr| — 6 04 
F5 = o-f - 10 af 0-2 + 20 af + 15 en al 
-30 CTi 0-4 - 20 02 0-3 -I- 24 CT5 



(7) 



and so on, in which we have set 



^ |0,|2™ = ^"-i/M. (8) 



The normalization of 0^ implies ai ~ 1. For the bulk exciton ground state, /(m) — 16 (Sm — 5)!!/(8to — 2)!!. By 
carefully counting the number of occurences of cycles of length 2 m leading to Um , it is possible to show that the 
general expression of F^v is 



F.=m^ n- ^ 

where the [im] are determined by ^ mim = N . An even more compact expression of F^ can be found 



(10) 

i)=0 



S(r/)=n(l + '7l<^^l')- (11) 

fe 

However, these expressions are of limited practical usefulness for A^ ^ 1: eq. (||) already contains more than 10^ 
different terms for A^ = 60, and S is a complicated function of 77. For numerical purposes, the recursion formula 

FN^f:('^r'^'w^,--p--m (12) 

m— 1 ^ ^ 

with Fq = 1 turns out to be more convenient. 
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We first note that eq. (0) implies > while, from eqs. Fn is a polynom of order iV — 1 in z. Since 

0-2 = 33z/128 = (33 7r/2)</T/, we find that F2 > for V^^^ = L > L2 ~ 3.7 a^:. For large N, numerical calculations 
show that Fn > for L > Ln ^ ia^ N^^^. As expected, excitons should be a few Bohr radii apart to be considered 
as bulk excitons. For smaller volumes, excitons distort and the modified accordingly. 

From eqs. (||)-(|^), we find that the z-expansion of Fn starts as 



F, 



N 



l-(72iV(7V-l)/2 + 0((iV2z)2), 



(13) 



leading to a N'^ dependence of the correction to 1. The exact value of a2 given above shows that Fn — 1 for 
25 N'^ a^/V ^1. A simple way to understand that this N'^ dependence comes from the iV-body character of Pauli 
exclusion is to return to the definition of . We have 



ki,k' 



X 



(14) 



The matrix elements in eq. (|lj) give ±1 (their product being +1) for {ki}—{k'j} and the ki all different; if not, they 
give 0. This allows to reexpress Fn as 



E 



felV^{fel.fc2.--- 



(15) 



1} 



As expected from Pauli exclusion, eq. (|T^) merely states that the A^*'' exciton must be built with a wave vector 
different from those already used by previous excitons. Without this restriction, all the sums would be equal to 1, 
and so would Fn, as for perfect bosons. 



The N'^ dependence of Fn may be obtained easily from eq. (15) using a model calculation, which contains all the 
physics of the problem. An exciton is basically an object of spatial extension a^, so that in a volume V = {L ^ Gx) 
the wave function in fc-space has a width l/a^, all the fc's being 2 tt/L apart (see Fig. 0(b)). Instead of its exact value, 
let us take (j)k — 4>q up to k — l/a^, and otherwise (see Fig. |^(c)). The normalization condition implies A^'o 4>q = 1, 
where Nq — [{1 / ax) / {2 tt / L)]^^ is the number of k's such that (pk ^ 0. For N < Nq we get 



Fn+1 = 1 (1 - ^l) 



{1-N4,l) 



(16) 
(17) 



For A^^ <C 1, eq. (p^ gives Fn ~ 1 - A^ (A^ - 1) (f>l/2 in agreement with eq. (Oh, whereas eq. (p7|) gives 



Fn 



^N{N~l) 4,1/2 



for < 1. 



(18) 



Equation (jlq) shows that although the excitons are "close" to bulk excitons for N (j)^ <^ 1, their normalization factor 
(0|i?^ {B'^)"\0) is much smaller than A^! (as for perfect bosons) if A^ < 1 < Af ^ 0o- This A^^ dependence may be 
viewed as a "stimulated" many-body effect. Indeed, Pauli exclusion being intrinsically a N-body interaction, many- 
body effects based on this interaction are quadratic. Note that this exponentially small factor Fn turns out to have 
a relatively small influence on A^! because 



If we now consider the exact 



NIFn eNilnN-l-{N-l)4'„/2+-)^ 

bk, we can actually show that Fn has a similar behavior 



F, 



N 



^-N [N-l] 0-2/2 



for AfCT2 < 1, 



(19) 



(20) 



in agreement with the first terms of the expansion given in eq. (^3|) . 

We now turn to the expectation value of D taken on l^**^^)}. It is convenient to introduce defined by [D, fit] 
2 Ct. Since [C^, St] = and DjO) 0, we get 



D{B^)'^ |0) = 2N{B'^)'^-^C^ |0). 



(21) 
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This leads to 



B (Bt)JV+i |o) = {N + l) (St)^ |0) 

-TV (iV + 1) (5^)^-1(7^10). 



From eqs. (^ij) and (|2^), we easily deduce 



(0|g^£)(gt)A^|0) 
(0|B^ (St)^|0) 



2 1 



Af+l 



F, 



N 



(22) 



(23) 



We check that {D)^ — for perfect bosons. For our model 0fc, eq. ( [iq ) immediately gives (-D)^ = 2N(/)q, while 
for the exact exciton wave function, we find 



2Na2 



for N(T2 ^l. 



(24) 



Taking (i?)^ <C 1 as the criterion for N bulk excitons to behave as bosons leads therefore to a prefactor A = 33 tt ~ 100 



Another relevant operator is the boson number in the |5'(^)) state. From eqs. (|l|) and (^ 
show that 

{B^B),=N-{N-1){D)J2. 



it is straightforward to 



Obviously, (B^ B)j^ — N for perfect bosons because (D), 
from eq. (|24| ) we get 



0. Its value is N — N {N — 1) (jj^ with our model 0fc 



(StB}„ - iV(l - (iV- 1)^2) forAr(72<l. 



(25) 
while 



(26) 



The percentage of pairs not counted as bosons is thus small for A^cr2 — 50 N a^/V 1. 

The above formalism can be readily applied to 2D excitons. From the 2D ground state wave function we have 
= 2//(l + a?/4)^, where y — 27r a^/S*, S being the sample's area and still the 3D exciton Bohr radius. We 
find am = 2y/(Zm— 1), so that a2 = 2y/5. Numerical calculations give L^v — l.ha,^ N^/'^ for N ■» I, while the 
condition 27V (T2 *C 1 now reads bNa^/S ^ 1. Assigning to the 2D exciton an effective radius a2D — ax/4, which 
corresponds to the maximum of the radial charge density [lO| ], we find again that our new criterion differs from the 
geometrical one N o^d/S <IC 1 by nearly two orders of magnitude. 

Let us end by reconsidering the possibility of Bose-Einstein condensation of 3D excitons in the light of these 
new criteria. N noninteracting bosons undergo Bose-Einstein condensation when their density n — N/V exceeds a 
temperature-dependent critical value nc{T), which, when applied to excitons, is given by H 



ksT 
R 



6.62 



TOe rrih ( al ndT) 



{me + muY 



9 



2/3 



(27) 



where R and g are the exciton binding energy and degeneracy, respectively. For CU2O orthoexcitons pd| |, R ~ 
150 meV, Ux — 0.7 nm, g = 3, and rrie ~ 1, mt ~ 0.6, so that lOOnca^ — 0.1 for a carrier temperature T — lOK, 
which means that for n > Uc, excitons barely obey our criterion 100 n aj^ ^ 1 for bosonic behavior. Besides, considering 
a typical excitation volume V ~ (10 //m)^ x (l/zm) necessary to get a large enough density, we find 25 V ^ 10^. 
On account of our criterion 25 A''^ a^/V <gC 1, Fn would be extremely small. It thus seems unlikely that excitons could 
be considered as noninteracting bosons under these experimental conditions. The critical density should therefore 
differ from eq. (^^. In contrast, let us stress that both our criteria are fulfilled in the case of atoms when we put 
numbers corresponding to actually observed Bose-Einstein condensation. Indeed, for qq ^ 0.1 nm and condensates 
of up to 5 X 10^ atoms, we find at most 100 N a^/V ~ 10~* and 25 A^^Oq/T^ ~ 10"'^ for the different configurations 
reported in the literature fl^ ]. The main reason for such a difference lies with the temperature T, which is about 
10~^K in atomic condensates. Excitons being "transient" particles produced by excitations, they cannot be cooled 
down as much. 



Conclusion 

While previous criteria for bosonic behavior of excitons were related to the Coulomb interaction between electron- 
hole pairs, we have studied the consequences of the fermionic character of these excitons induced by Pauli exclusion 
between their constituents. They lead to prefactors two orders of magnitude larger than the usual Mott criterion 
N a^/V ~ 1, so that excitons can be considered as bosons up to densities well below the Mott density only. This is of 
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importance for the possible observation of the Bose condensation of excitons. We have also found that (0|i3^ {B^)'^\0) 
is exponentially smaller than its boson value VV! unless 25 N"^ a^/V <C 1. This unexpected iV^ dependence has to be 
related to the intrinsic A''-body character of Pauli exclusion. 
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FIG. 1. Representation of the ground-state excitonic wave function in real space (a) and fc-space (b). The model wave 
function is represented in (c). 
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Figure 1 
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